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Abstract 

After a brief discussion of the Bogoliubov inequality and possible 
generalizations thereof, we present a complete review of results con- 
cerning the Mermin- Wagner theorem for various many-body systems, 
geometries and order parameters. We extend the method to cover 
magnetic phase transitions in the Periodic Anderson Model as well as 
certain superconducting pairing mechanisms for Hubbard films. The 
relevance of the Mermin- Wagner theorem to approximations in many- 
body physics is discussed on a conceptual level. 

1 Introduction 

The quest for criteria for the existence, or absence, of phase transitions 
in physical systems has been a dominant theme in theoretical physics ever 
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since the phenomenological concept of a phase was introduced in the con- 
text of equilibrium thermodynamics which is associated with the names of 
Bohzmann, Gibbs, Ehrenfest and many others. The pioneering work in the 
phenomenological theory of phase transitions was soon supplemented by at- 
tempts to provide a detailed theoretical description of macroscopic phenom- 
ena via microscopic many-body theories. The Lenz-Ising model (1925) was 
the first such model, and Onsager's solution of the 2-D Ising model (1944) 
is evidence that, by then, the mathematical physics of phase transitions had 
become a subject in its own right. Other landmark developments of the 
1950s and 1960s were the famous Yang and Lee papers (1952) providing a 
mathematically rigorous scenario of how phase transitions can occur in the 
thermodynamic limit, and the papers by Hohenberg (1967) and Mermin and 
Wagner (1966) excluding certain phase transition in low-dimensional systems. 

It is the latter result, generally known as the Mermin- Wagner theorem, 
which we shall discuss in this paper. In doing so, our goal is three-fold. First, 
we discuss apphcations, and possible generalizations, of the Bogohubov in- 
equality, which underlies the proof of the Mermin- Wagner theorem. Then we 
present a survey of existing proofs of the absence of finite-temperature phase 
transitions in low-dimensional systems along with recent generalizations. Fi- 
nally, by means of a simple example, we discuss the relevance and validity of 
the Mermin- Wagner theorem in approximate theories designed to describe a 
magnetic phase transition. 

2 Mathematical tools 

The occurrence of a phase transition is often intimately related to the failure 
of one of the phases to exhibit a certain symmetry property of the under- 
lying Hamiltonians. Crystals, for example, by their very lattice structure, 
break the translational symmetry encountered in the continuum description 
of fluids; ferromagnets, in addition to the spatial symmetry-breaking due to 
their crystal structure, are not invariant under rotations in spin-space, even 
though the underlying Hamiltonians describing the system may well be. Less 
obvious types of symmetry-breaking occur in other quantum systems, such 
as superfluids and superconductors, where a breaking of gauge invariance 
occurs. 
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2.1 Bogoliubov quasi-averages 

Bogoliubov has devised a method for describing the occurrence of spon- 
taneous symmetry-breaking in terms of gwasz-aferag'es. Normally, for 
systems in statistical equilibrium, the expectation value of an operator A 
is given by the trace over the equilibrium density (or statistical) opera- 
tor p = exp(— /97i) times the observable A. Thus, for the infinite system 
(y ^ oo) one must calculate 



(A) = lim tr{pA) 

V^oo 

lim tr [e-f^^A) , (1) 



where 7i is the grand-canonical Hamiltonian, 7i = H — pN (N: number 
operator). However, it turns out that under certain conditions such averages 
may be unstable with respect to an infinitesimal perturbation of the Hamil- 
tonian. If a corresponding additive contribution H^, = uH' of the order of v 
is added (where u is a small positive number which will eventually be taken 
to zero: z/ ^ 0), i.e. 

H, = H + H,- pN, (2) 
one can define the quasi-average of A in the following way: 

(A)^ = lim lim tr^e''^'^" A). (3) 

The average ([^) and the quasi-average (0) need not coincide, since the two 
limits in (|^) may fail to commute within some parameter region (i.e. for 
some combination of p and (3). 

Quasi-averages are appropriate for cases in which spontaneous symmetry- 
breaking occurs, as can be shown by a simple argument. Suppose the Hamil- 
tonian Ti displays a continuous symmetry S, i.e. it commutes with the 
generators Tg of the corresponding symmetry group, 

[^,ry_ = o. (4) 

If some operator B is not invariant under the transformations of S, 

[s,ry_ = cVo, (5) 
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the (normal) average of the commutator C* vanishes, 



(C) = 0, (6) 

as can be readily seen from eqn. ([^) by use of cyclic invariance of the trace. 
In those instances, however, where the perturbative part H^, of (|]) does not 
commute with F^, this will give a non-vanishing quasi- average: 

= limtr(e-^«^[5,ry_) 

^ 0. (7) 

Thus, even though one might naively expect the quasi-average to coincide 
with the "normal" average in the limit u ^ 0, quite generally this will 
not be the case. Note that the quasi-average depends on the nature of the 
perturbation added to the "original" Hamiltonian. 

As a first example, let us examine the Heisenberg model 

which is invariant under the continuous rotation group generated by the to- 
tal spin vector S = ^ ■ Si because of 



H,S 



0. Thus, we can take S as 



the operator B, and from (j^) one finds [S", S''^] _y = 0, where S'",S''^ are 

the components of the total spin vector S. Together with the commutation 
relations for spin operators, e.g. [S^, S^]_ = ihS^, it is obvious that the con- 
ventional average of the magnetization vanishes. This is just a manifestation 
of the fact, that on the macroscopic level, for an ideal, infinitely extended sys- 
tem, there is no preferred direction in space. One can think of this situation 
as a "degeneracy" with respect to spatial orientation. Adding an external 
field. Bo = BqCz along the z-axis for example, lifts this degeneracy and via 
the contribution Hi, ~ BqM{T, Bq) to the Hamiltonian (M: magnetization) 
one can, thus, construct appropriate quasi-averages for arbitrary operators 
according to equation (H). 



2.2 Order parameters 

In the theory of phase transitions, the first step is to identify a quantity 
whose (quasi-) average vanishes on one side of the transition, but takes on 
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a finite value on the other side. This quantity is called the order parame- 
ter. In a continuous phase transition, the order parameter may gradually 
evolve from zero at the critical point to a finite value on one side (usually 
the low-temperature side) of the transition. For different kinds of phases, 
different order parameters must be chosen. From a phenomenological point 
of view, one must consider each physical system anew. We have already 
mentioned the liquid-vapour and the ferromagnetic-paramagnetic transition 
as prototypes for phase transitions. In the former, the obvious choice for 
the order parameter would be the difference between the mean densities, i.e. 
P — Pvapour', in the latter the relevant order parameter is the magnetization 
M. Within a given many-body model, the magnetization can be defined in 
microscopic terms, as we shall see shortly. 

Sometimes, as in the transition to the superconducting state, it may even 
be possible to characterize the same type of phase transition by use of dif- 
ferent order parameters. According to the standard theories of superconduc- 
tivity, at low temperatures electrons with opposite spins form Cooper pairs. 
Thus, a possible order parameter would be the average probability amplitude 
to find a Cooper pair at a given lattice site in the crystal. Alternatively, one 
could characterize the phase transition through the gap parameter whose 
modulus is the difference in energy per electron of the Cooper pair con- 
densate and the energy at the Fermi level. We shall briefly return to the 
problem of superconductivity when we consider the possibility of pairing at 
flnite temperatures in fllm systems. 

2.3 Bogoliubov inequality 

The Bogoliubov inequality is a rigorous relation between two essentially ar- 
bitrary operators A and B and a valid Hamiltonian H of a physical system. 
In its original form, proposed in [|l|, it is given by 

\{[C,Al)\' <^{[A,A^]^){[C\[H,CU_), (9) 

where (3 = l/ksT is the inverse temperature and (....) denotes the thermo- 
dynamic expectation value. A and B do not necessarily have an obvious 
physical interpretation from the very beginning, so the physical signiflcance 
of (P) will depend on the suitable choice for the operators involved. The 
inequality can be proved by introducing a scalar product which is based on 
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the energy eigenvalues En and the orthogonahty of the corresponding en- 
ergy eigenstates |n) of the Hamiltonian H and to which one then apphes the 
Schwarz inequahty. The details of the proof are readily accessible elsewhere 
(see e.g. textbook [^); however, we briefly note that as a result of this deriva- 
tion, the two factors on the RHS of (P) each are upper bounds to a norm and 
can, thus, be bounded from below by zero. In particular, if, for example, the 
double commutator depends on some parameter k, one will always find 



[iC HU C^]) (k) + ( [IC HU C']) {k') > ( [IC HU C^]) (k). (10) 



Dividing both sides of (P) by the double commutator and summing over 
all wave vectors k associated with the first Brillouin zone in the reciprocal 
lattice, one arrives at 



[[[C,H].,Cn,){k) 



(11) 



2.4 Generalized Bogoliubov-type inequalities 

The use of inequalities in the theory of phase transitions has developed into a 
subdiscipline of its own right within the field of mathematical physics 
consequently, this article is not intended to give a complete account of such 
approaches. The close kinship between some of these methods and the Bo- 
goliubov inequality, however, justifies a brief discussion of what one might 
call generalized Bogoliubov-type inequalities. The common feature of these 
inequalities is the fact that certain algebraic properties of the quantities in- 
volved allow one to use the Cauchy-Schwarz inequality, thus leading to upper, 
or lower, bounds for physical observables, such as the dynamical structure 
factor 



CBAik',k; E') 



^^^^E't/n J2 in e*^*/^5(^', Q)e-''^"''A{k, 0) 



n 

NE 



m \A{k,0) 



n 



n){n 



(12) 



5(fc',0) m)5{E' -{En,-En)) 



6 



(where the sum runs over all energy eigenstates |n) , |m) corresponding to 
eigenvalues En, Em), or the susceptibility 

OO 

1 f ,^,{^-^~ncBA{E') 



XAB{k,k';E,ir]) = — dE 



2tx 



E-E' + ir] 



Wagner H discusses the static susceptibility Xab, defined as 



XAB{k,k') = yi[ Urn XAB{k,k';E,ir]) 

E,ri~*0 



V / dE 



,SA B{k,k'-E') 
E' 



(13) 



in terms of the spectral density which, in its spectral representation, is given 
by 



h 



SAB{k,k'-,E') = ^Y.kA'^^^^^) 



m } { m 



B{k', 0) I n) e-^^- (e'^^' - l) 5iE' - (E^ - E^)). 

nm 

(14) 

It can then be shown 0,14^ that 

'A{k)-B{k')) :=XAB{k,k') (15) 

is a valid scalar product in the space of operators A, B. The Schwarz in- 
equality for the static susceptibility then reads 



XAB{k,k') <XAAAk)XBw{k')- 



(16) 



Physically significant relations can be deduced from this general statement 
by an appropriate choice of operators A and B. If A is taken to be a time 
derivative of another operator H, i.e. 



the inequality 



V I dE^^^^^ > 
E 



(17) 



(18) 
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follows, which relates the response function for the observable B to commu- 
tators that can, in principle, be calculated directly from one's knowledge of 
Q, B, and H. If for a given many-body Hamiltonian H one specifies Q and 
B further, it is possible to obtain bounds for non-trivial order parameters. 
Thus, for the planar magnetization Mp defined as 

Mp(T, Bo) := ^ E + (19) 

i 

(where i runs over all lattice sites, erf and erf denote the respective compo- 
nents of the spin, and the real constants and f3i can be chosen so as to take 
into account various kinds of spin ordering within the xy-plane) , evaluating 
([I8| ) with operators Q and B given by 



and 

j 

gives an upper bound for Mp in terms of the longitudinal susceptibility x^^{k)- 
Since for the latter, rigorous limits and estimates have been established, e.g. 
for the attractive Hubbard model 0, these carry over to the order parameter 
Mp. 

Before turning to exact results concerning finite-temperature phase tran- 
sitions, let us discuss how correlation inequalities can be put to use in the 
zero-temperature case. Pitaevskii and Stringari have suggested to define 
a scalar product simply through the anticommutator: 

{A; B) := {[A\ BU) (20) 

which gives the Schwarz inequality 

{[A\AU){[B\BU)>\{[A\BU)\\ (21) 

It is then possible to define auxiliary operators, denoted by a tilde, such that 

(n C 0^ := {n\C\0) 

(o C ny.= - {^\C\n) (22) 
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from which it follows that 



n 







(7t 



n 



^ = _<^|ct|0) 
= {Q\C^\n). 



(23) 



At T = the expectation values in (^T]) are evaluated in the ground state 
|0) , and if instead of B the newly defined operator B is used, one deduces 
from 



[A\B]^ 
[B\Bl 



(0|[At,i?]_|0> 
(0|[5t,5] + |0> 



(24) 



that the Schwarz inequality holds not only for the anticommutator on the 
RHS but also for the commutator, so that (writing down ( [2T| ) for the opera- 
tors v4, 5, and making use of (p^ to express B in terms of B) one finds 

|2 



'\A\A]^){[B\B]^)>\{\A\B]^)\\ 



(25) 



With this in mind, one can turn to the dynamical structure factor, which for 
the transverse spin-spin correlation function at T = is given by 



C+-(g,E) 



T=0 



1 



lT=0 



cite^^*/^(a+(g,t)a-(-g,0)>| 
^ j t/te*^*/'^5^(0|e*«*/V+(g,0)e-*^*/'^| 



-oo 

CO 



'X 



J2 - (^,n - ^o)) I (0 |a+(g, 0) \m)\\ 



Since Eq is the ground state energy, we have Em > -Z?o, so the ^-function only 
gives a contribution if > 0. As a corollary 



EC^-{q,E) 



> E > 
= ^ E < 0. 



(27) 



If one defines 



C^{q,E) = -{C+-(q,E)+C-~'i-q,E)) 
= ^([a+(g,E),a-(-g,E)]_ 



(28) 
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and 

C^{q) = j dEC^{q,E), (29) 

one can make use of P3| ) for operators cr+(g + Q), a^{—q) to get 

C^{q + Q)C±{q) > \m{Q)\\ (30) 
where the magnetization 

M{Q) = ^-^ {Y.<yte-'^'^^ (31) 

has been introduced (the reciprocal lattice vector Q as usual accounts for 
possible antiferromagnetic order). Eqn. (^Tf ) tells us that C±{q,E) is nonva- 
nishing only for E > 0, where it is positive, so the integral quantity C±{q) 
may be formally written as 



C^iq) = I dE^EC^E) . (32) 

In this form, one can apply Holder's inequality which states that for real- 
valued function / and g for which |/(x)|^ and |5'(x)|'^ are integrable (where 
p and q are numbers satisfying the relations + q^^ = 1 and p > 1) the 
following inequality holds: 

''dxf{x)g{x) < ( [\x\f{x)A ( f dx\g{x)A . (33) 



For p = q = 2 and with /, g defined as 

fix) = ^ECME) 



9i^) = \l ^ (34) 



Holder's inequality (^) applied to (^) gives 



C^iq) <Jj dEECM E).Jj dE^^^. (35) 
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The first of tlie integrals can in many cases be bounded from above by 



dEEC±{q,E) 



a+iq,E),[n,a~i~q,E)]_ 
f{q)Nq^ < const ■ Nq^ 



(36) 



which suggests writing 



/ dEEC^{q,E) 
JdE 



C_L(g,g) 

E 



= ^^(^) ■ 



thus defining the new quantity uj^q). Inequality 
using (|35D and rearranged by use of (|37D to give 



C±{q + Q) > 



M{Q) 



(37) 

can be strengthened 
(38) 



This equation relates the order parameter M, the transverse susceptibihty 
X~^~ and the structure factor (via C±). The relation between these 

quantities can be made even more pronounced by summing both sides of 
(|38| ) over the first Brillouin zone. One can then go on to show that in the 
thermodynamic limit a relation of the form 



(2vr) 



d'^q 



M{Q) 



l.B.Z. 



(39) 



holds (where turns out to be a model-dependent constant). Provided 
that uj{q) and X~^~{(l) remain finite in the limit g — 0, this means that in 
ID not even at zero temperature can there be a transition to a state with 
nonvanishing M{Q), as this would lead to a logarithmic divergence of the 
RHS and thus would contradict being a constant. 

It seems to us that the range of applications of correlation inequalities 
has not yet been fully exhausted. Using the mathematical tools discussed in 
the previous paragraphs, and variations thereof, should serve as a starting 
point for the derivation of similar relations for other order parameters and 
many-body models. 
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3 The Mermin- Wagner theorem: A survey 



Having established the wide range of relationships the Bogoliubov, and re- 
lated, inequalities can generate, we shall now turn to a narrower range of ap- 
plications, namely results concerning the absence of finite-temperature phase 
transitions in low-dimensional systems. 

Hohenberg was the first to note that the Bogoliubov inequality could 
be used to exclude phase transitions, showing that there could be no finite- 
temperature phase transition in one- and two-dimensional superfiuid systems. 
At roughly the same time Mermin and Wagner, following a suggestion by Ho- 
henberg, considered the case of spontaneous magnetization in the Heisenberg 



model. |T0[ Since Mermin and Wagner's proof has become the exemplar for 
studies concerning the absence of phase transitions, we shall briefiy outline 
the paradigmatic procedure (see also fill). The general idea is to use the 
Bogoliubov inequality (|ll|) to find an upper bound f{BQ,M) for the order 
parameter in question, e.g. the spontaneous magnetization M : 

M<f{Bo,M). (40) 

As indicated, the upper bound will normally depend on the external (e.g. 
magnetic) field that couples to the order parameter, and (implicitly) on the 
order parameter itself. To answer the question whether or not a phase tran- 
sition to a state with a non-zero value of the order parameter occurs, one 
must consider the case Bq 0, i.e. the behaviour of the upper bound in the 
case of vanishing external field. The subsequent argument against a phase 
transition proceeds by reductio ad ahsurdum: If the assumption M 7^ can 
be shown to lead to a violation of ( ^OD in the limit i?o — >■ 0, where equation 
(|40| ) is derived from the Bogoliubov equation (which is known to hold for the 
corresponding many-body system), it must be dropped. This, then, leaves as 
the only conclusion that Mjj^^O, the case of vanishing order parameter and 
no phase transition. This argument, of course, only succeeds if the initial 
inequality is indeed true "a priori, " and the most straightforward way to 
achieve this is to resort to the Bogoliubov inequality. From this it follows 
that, once one has specified a many-body model by its Hamiltonian if, the 
operators A and C in the Bogoliubov inequality must be carefully chosen so 
as to give the desired order parameter. 
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As an example, consider the Heisenberg model 

H^He^) = - Yl J^AS^Sr + S^S])-b e^'^'^'S,^ (41) 

ij=l...N 2=1. ..Af 

(where the term b'^-e'^^'^'S^ is due to the interaction with an external 
magnetic field b = xhe relevant order parameter is the magnetiza- 

tion 

M = ^'-^T.'-'''''' {SD (42) 

i 

where the factor e~*^'^* already accounts for antiferromagnetic order by 
changing the sign of the spins in one sublattice (given that K has been 
properly chosen to achieve just that). This suggests that an upper bound for 
{S^) is essential. Comparison with ([TTD suggests the choice 



A = S~{-k-K) (43) 
C = S+{k) (44) 

as, by virtue of the commutation relations, 

's''ih),S'ih)]_ = 2hS:^{h + h), (45) 

this will indeed give a contribution proportional to the magnetization on the 
left-hand side of the Bogoliubov inequality: 

{[C,A]_)=^,NM{T,Bo). (46) 

(^(j) denote constants depending, at most, on fixed parameters of the many- 
body model, i.e. in this case the Heisenberg exchange integrals.) The other 
(anti-)commutators that feature in the Bogoliubov inequality can be bounded 
from above, as has been shown by Mermin and Wagner: 

[[C,H]_,C^_) < e2N{\BoM{T,Bo)\+^sP) (47) 
J2([A,A^]\ < 2UN\ (48) 
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For the Bogoliubov inequality, 



k k 

one finds in this case 

I 55iv(|BoM(r,B„)| + 6Pj 

In the thermodynamic hmit, the only situation where one can at all hope 
for a phase transition, the sum is to be replaced by an integral, e.g. for the 
two-dimensional case 

k 

where L^/(27r)^ is the area in fc-space associated with one quantum state. 
Restricting the support of the integral to a finite-volume sphere inscribed 
into the first Brillouin zone only strengthens the inequality, so 

fk] L^M\T,Bo) f dk- — ^- -<e4/3 (51) 



where /cq is the cutoff corresponding to the sphere in /c-space. In the thermo- 
dynamic limit under consideration, L and approach infinity in such a way 
that the specific volume v^'^ = L?' /N remains finite throughout. Evaluating 
the integral, one then arrives at 

M\T,B,)<i— \ (52) 

V ^ \BqM{T,Bo)\ ) 

As -Bq ^ 0, the denominator diverges logarithmically, thus forcing the mag- 
netization to vanish. This result is independent of the original choice of the 
auxiliary wave vector K (see (^2])), so that both ferromagnetic and anti- 
ferromagnetic order is ruled out in the two-dimensional Heisenberg model. 
One easily verifies that a similar divergence of the denominator rules out 
spontaneous magnetic order in the one-dimensional case. 
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The main steps in any proof of the Mermin- Wagner type are: a) choice of a 
many-body model characterized by its Hamiltonian if; b) identification of the 
order parameter for the phase transition to be discussed; c) adequate choice of 
operators A and B in the Bogohubov inequahty so as to single out the order 
parameter; d) search for non-trivial upper bounds for the (anti-) commutators 
in the Bogoliubov inequality; e) proof that in the thermodynamic limit the 
assumption of a spontaneous non-zero value for the order parameter will be 
self-refuting in the one- or two-dimensional case (and, possibly, other cases as 
well). This basic scheme has been applied to a wide range of different models 
and order parameters. The remainder of this section attempts to present a 
survey of established as well as new results, along with brief discussions and 
related references. 

Shortly after the papers by Hohenberg |Q and Mermin and Wagner flO 



Mermin showed that a classical Bogoliubov-type inequality holds, if in equa- 
tion d^) one replaces the commutators by Poisson brackets and also requires 
certain quantities in the corresponding classical thermal averages to vanish. 
For short-range interactions, one can then rule out phase transitions for the 
"classical," i.e. infinite-spin limit. []T2| 



Walker used the Bogoliubov inequality to rule out the possibility of a 



phase transition, in one and two dimensions, to an excitonic insulating state. [13 
The main idea behind the formation of the insulating state involved the as- 
sumption that certain conditions might favour bound electron-hole pairs to 
form if a semiconductor with a very small gap, or a semimetal with a very 
small band overlap, was to be cooled to a sufficiently low temperature. The 
proof that such a state is impossible in the low- dimensional case is straight- 
forward for electrons interacting via a potential which falls off faster than 
\f\^^ {D: number of dimensions); it can be generalized to cover certain sim- 
plified Hamiltonians with weaker convergence behaviour. A generalization 



to an isotropic two-band model can also be given. ||13 



The Bogoliubov inequality, together with an analogous classical inequal- 
ity, also rules out the possibility of crystalline ordering in two dimensions, 
thus confirming earlier suggestions by Peierls and Landau that there could 
be no two-dimensional crystalline ordered state. [1^1 This is significant, since 
the extension to crystalline ordering is not quite straightforward: Contrary 
to types of ordering where the energies of fiuctuations that cause the disorder 
are kinetic (e.g. superfiuid, superconducting, or excitonic interactions), the 
relevant energy in the crystalline case is potential. Also, one cannot simply 
posit that one particle only interacts with a finite number of neighbours (as 
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one can in the case of spin systems; see next section), because a given particle 
may well diffuse through the crystal and interact with many, and possibly 
all, others. The proof can nevertheless be achieved by assuming a system of 
identical particles with pair potential $(r), where $ as well as the related 
function 

^(f) =$(f)-A|f|^|A<l>(f)| (53) 
are required to satisfy the criteria for the existence of a proper thermody- 



namic limit for a sufficiently small positive value of A.|15] While this includes 



Lennard- Jones type potentials, it does not cover hard-core potentials. 



m 



The corresponding result for the quantum case, only sketched in |n 
an appendix, has been reproduced in detail by Fernandez. |]16| The paper dis- 



cusses an electrically neutral sytem of nuclei and electrons which is confined 
to one or two finite dimension(s) with the remaining two or one dimension(s) 
being infinite. The paradigmatic examples for such geometrically restricted 
systems are the infinite slab and the rod with rectangular cross section and 
infinite length (see section 3.2 below for further examples of this kind). One 
then proceeds to show that for systems of this kind, no maximum long-range 
crystalline order can exist, or, formally: 

(^^.) = for any K 0, (54) 

where 

i^g) = hm (p^) /N (55) 

with Pg = J dfp{r) exp (^—iK ■ (56) 

and p(r) is the number-density operator, and the limit is taken with N/V 
being held constant (V^: volume of the system). 

Kishore and Sherrington considered a quite general Hamiltonian of 



electrons interacting non-relativistically among themselves, and with spa- 
tially ordered or disordered scatterers, excluding spontaneous low-dimensional 
magnetic order. The restriction to non-relativistic interactions means that 
spin-orbit effects are excluded from the problem, thus not ruling out a phase 
transition for Ising-type, or other suitable anisotropics. The role of disorder 
and impurities has also been discussed by Schuster, |T^] who discusses the 
influence of a static random field conjugate to the order parameter. Such 
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random fields might be generated by electrically charged or magnetic impu- 
rities in certain quenched systems. The Bogoliubov inequality, in that case, 
allows one to exclude magnetic order for a classical X — Y model with Gaus- 
sian distributed random field in less than or equal to four dimensions. The 
proof makes use of Mermin's classical analogue of the Bogoliubov inequality 



mentioned earlier in this section. [12| The significance of the Mermin- Wagner 



theorem for the X — Y model and its characteristic Kosterlitz-Thouless- 



Berezinskii transition is discussed in 19 . For the X — Y model it has also 



been pointed out |]20[ that finite-size magnetization may obscure Mermin- 



Wagner-type behaviour for any realizable physical system. 
3.1 Spin systems 

Much research has been carried out in order to extend the statement to spin 
systems other than the simple Heisenberg model considered by Mermin and 
Wagner. Wegner |2T[] considered a model describing a system with locally 
interacting itinerant electrons with pseudospins a, to which the i?o-field cou- 
ples, giving a contribution ~ Bq ( 1 to the Hamiltonian. 
The dynamics of the electrons, in this model, are governed solely by the ki- 
netic energy Te = ^gg. (^ka- '^^^ interaction V of the electrons with the 
nuclei and among themselves, is further assumed to satisfy the condition 

[C,V]^ = (57) 

(where C is the operator that appears in the Bogoliubov inequality). 

A more realistic model has been discussed by Walker and Ruijgrok [Q: It 
includes Coulomb and exchange effects, with possible non-local interaction. 
The authors consider a model for which Lieb and Mattis had previously ruled 



out ferromagnetic ordering in one dimension EH] , and which turns out to be 



a special case of their general many-band model for interacting electrons in a 



metal. Ghosh ||2J], more specifically, recovered the Mermin- Wagner theorem 



for the Hubbard model given by 

Huub = X] X] Tijcj^Cj^ + ^ 5Z '^i'^^i-'^ + ^0 X](^iT - ^u) exp(-zi? ■ Ri). 

ij cr i I 

(58) 

Van den Bergh and Vertogen present a proof of the Mermin- Wagner 
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model for the s-d interaction model characterized by the Hamiltonian 



Ima Im 



(59) 



Just as in the case discussed by Walker and Ruijgrok, the s-d interaction 
does not satisfy Wegner's criterion (|57D . Furthermore, the system is com- 
posed of two subsystems, so that one has to show that both the conduction 
electron system and the system of localized magnetic moments fail to exibit 
a magnetic phase transition in ID and 2D. 

Robaszkiewicz and Micnas have generalized the proof for the s-d model by 
including interactions within the spins of each subsystem and Hubbard-type 
interactions; the model, thus, covers the modified Zener model, the extended 



Hubbard model and s-d models as particular cases. p6| 

The range of validity of the Mermin- Wagner model can, of course, also be 
extended by discussing more general geometries: Baryakhtar and Yablonskii 
P7| , for example, have shown that the Mermin- Wagner theorem remains valid 
for systems with an arbitrary number of magnetic sublattices, and also ex- 
cludes non-collinear magnetic order when an external field is applied. Thorpe 



| 28| considers the case of ferromagnetism in phenomeno logical models with 



double and higher-order exchange terms; similar results were obtained in the 



multi-sublattice case by Krzemihski.[29| These methods differ from others 



in that the Hamiltonian is written as a series expansion in terms of spin 
spherical harmonics, thus offering a more systematic way of evaluating the 
Bogoliubov inequality using the defining properties of spin spherical har- 
monics. A closely related proof, using spherical tensor operators, has been 
put forward for the problem of ordering in quadrupolar systems of restricted 
dimensionality. |]3D[] 

Uhrig has shown that, at finite temperatures, there can be no pla- 
nar magnetic order in the one- and two-dimensional generalized multiband 
Hubbard model 

H gen. Hub ^ ^ J^ici,j"/Cj^aa^jy<^ ~^ ^ ^ UiaJyTT'ia^jy (^0) 

ija ij 

ay ay 

Planar magnetic order, in this case, is taken to be characterized by an order 
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parameter of the form 



ia 

where crj^'^'^"* denote the pseudo-spin components at lattice site z, a is the 
band index, and rjia, Qa are essentially arbitrary real constants that fix the 
direction and the norm of the order parameter Mpian- Except for the choice 
of operators C and A in the Bogoliubov inequality, the proof of the Mermin- 
Wagner theorem carries through in the usual way. 

Attempts to generalize the Mermin- Wagner theorem to anisotropic mod- 
els with n-th nearest neighbour exchange interactions produce results that 
seem somewhat artificial. [^] Similarly, attempts to extend the Mermin- 
Wagner theorem to special cases in three dimensions remain somewhat incon- 
clusive. [0311 In both certain circularity enters the argument: the pa- 



rameters of the respective models are chosen "ex post facto" so as to enforce 
the absence of a phase transition. 

A comparatively large class of many-body models has been discussed by 
Proetto and Lopez. 0] They confirm the Mermin- Wagner theorem for Ander- 
son and Kondo lattices, which is of interest since the nature of the exchange 
interaction e.g. in the Anderson model is quite different from the Heisenberg 
model and variants thereof. The effective exchange interactions are higher- 
order functions of the hybridization matrix elements and are mediated by 
the conduction band. As can be shown via a canonical transformation, at 
higher order multisite interactions exist that go beyond the simpler models 
discussed thus far. The full Hamiltonian for this case is given by 

HAnd = y ] ^ij^ia^j<T + ^/ ^ ^ fiafio- + ~^ ^ ] fiafio-fi-afi-o- + ~^ ^ ] ^L^«o-cJ_, 
ija i(T i(T i(T 

( E/-/- ) ) + E (^(^^ - + ^-^O 

i \ a / \ a' / ijcr 

~Y E ((■^/t'^'T - fiih) + (cltCiT - cli^n)) exp (^-iK ■ Ri^ , 
I 

thus covering Hubbard model, Anderson model, Falikov-Kimball model, and 
Kondo-lattice model as special cases. The Anderson model has also been dis- 
cussed by Noce and Cuoco who rederive the Mermin- Wagner theorem for the 
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magnetic phase transition. Their paper also presents analogous proofs for 
different pairing mechanisms and the corresponding superconducting phase 
transitions. For the Hubbard model, a proof of the absence of superconduct- 
ing long-range order for different kinds of pairing has been given in ref.|]36[ 



A generalization of these cases will be discussed in section 3.3 below, along 
with its implications for superconductivity in thin films. 

3.2 Partially restricted systems 

Early after the publication of the papers by Hohenberg, Mermin, and Wagner, 
some authors had expressed their doubt as to the applicability of the Mermin- 
Wagner scheme for sytems of finite cross section (or finite thickness). As 
we have seen in the previous section, the scheme based on Bogoliubov's 
inequality ultimately rests on the divergence of the integral 

k<ko 

in one or two dimensions {D = 1, 2). The doubts were founded on the obser- 
vation that the wave function of a system in a box of finite cross section must 
vanish on the walls, thus forcing the wave function to have some non-zero 
curvature even in the ground state. This is not compatible with the wave 
function assuming a uniform value - which, it was believed, should corre- 
spond to the k = value ( "vanishing momentum" ) that gives the divergence 
in the first place. Chester et al. have argued that these doubts are 



ill-founded, since the fc-value is a purely mathematical index which need not 
be identified with any physical property of the system - just as in (|l^) the 
physical meaning of the /c-values is irrelevant to the mathematical inequality. 
This observation indicated that it might well be possible to extend the scope 
of the Mermin- Wagner theorem to partially restricted systems as well. 

Costache and Nenciu |^8[ discuss the overall magnetization of a partially 



finite three-dimensional Heisenberg model and reproduce the generalization 
expected from the arguments by Chester et al. A similar, though much more 
exhaustive, discussion has been given by Fisher and Jasnow in two papers 



39[,^2l5 where the authors discuss Bose systems with respect to off-diagonal 
ordering, and spin systems with respect to a magnetic phase transition. The 
general procedure is to embed the space Q occupied by the physical system 
into an enclosing "box" A and, furthermore, to allow for a decomposition 
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of f2 into "sub domains" T. This makes it possible to distinguish different 
boundary conditions - such as a possible "corridor" surrounding Q. This is 
of interest for the discussion of explicit bounds on e.g. the spin-spin corre- 
lation function, but also on the limiting behaviour of e.g. surface contribu- 
tions in partially finite systems as the system approaches the thermodynamic 
limit. The "global" approach to partially restricted systems, i.e. the pro- 



cedure of slicing up physical space into the system's domain, its subdomains, 
and additional auxiliary "boxes", however, appears to obscure the more ba- 
sic question whether or not within a given many-body system (e.g. of film 
geometry) a phase transition can occur. This, as will be seen in the next 
section, is quite independent of which numerical values certain model pa- 
rameters display in certain subdomains (apart from some general symmetry 
requirements.) 

For the Heisenberg model, Corciovei and Costache have discussed 



the role of boundary conditions imposed on partially finite systems, focussing 
on boundary conditions of the "pinning case". The authors also discuss a 
classical analogue of the Bogoliubov inequality (which, it appears, is derived 
independently of Mermin's paper [l^), which allows an alternative proof of 
the non-existence of magnetization in the classical case. The case of the 
Hubbard model of thin films is mentioned by Sukiennicki and Wojtczak |^ 



unfortunately no explicit proof of the absence of a phase transition is given 
in their paper (contrary to what one of the authors suggests in ref. p4|). 

3.3 Systems with film geometries 

The problem of geometrical restriction imposed on physical systems, can be 
made much more explicit, by defining the many-body system as a film system 
from the start. Thus, no restrictions need to be imposed "by hindsight" , and 
any further specifications of interactions etc. can be implemented on the level 
of the many-body model itself. A proof of the Mermin- Wagner theorem for 



systems with film geometries has recently been given |]TT| for the Heisenberg, 
Hubbard, s-f and Kondo-lattice models. In this section, we briefly outline 
the procedure, thereby extending the validity of the Mermin- Wagner system 
to the Periodic Anderson Model with fllm geometry (where we shall rule out 
a finite value of the layer magnetization), and to a superconducting pairing 
mechanism in Hubbard films. 

The film geometry is incorporated into the many-body Hamiltonian by 
assigning each atom in the sample a double index (^,7), where n represents 
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the vector Rn of the underlying 2D-Bravais lattice and the Greek index 7 
identifies which layer of the film is considered. We shall assume that the film 
systems consists of d identical layers stacked on top of each other; the total 
number of lattice sites within one layer is A^. 

For the Periodic Anderson Model the film Hamiltonian then is as 
follows: 

(4./-. + fLc...) -bJ2e-'^-^^ + 4J(64) 

iacT ia 

where the c^"^^- and f^'^^ -operators denote fermionic annihilation (destruc- 
tion) operators for the electrons of the conduction band and the /-electrons, 
respectively. (The number density operators n{^^ = f-aafiaa and the z- 

components of the pseudo-spins, cr^/j-.^ = (^/2) (^n^l/^ — nfj^^ , are defined in 

the standard way.) 

For the Bogoliubov inequality the operators A and C are chosen as 

A7) = cj-{-k-K) + aJ^i-k-K) (65) 
^ = EC'/3 = E(<(^) + 4(^))' (66) 

(3 f3 

SO that A is layer-dependent, whereas C is layer-independent (the sum 
extends over the whole sample.) The (anti-)commutators appearing in the 



Bogoliubov inequality (11) can now be evaluated by making extensive use 
of the (anti-)commutation relations for the conduction electrons and the /- 
electrons. The Hamiltonian-independent quantities can be calculated in a 
straightforward way: 

) = M,iT,Bo) (67) 

where the layer-dependent magnetization has been introduced (again with 
a phase factor e~*^'^" to account for possible antiferromagnetic ordering): 

M,,(T, B„) = 5^ e-" '^' + . (68) 



22 



The right-hand side of ( ]TT| ) can be bounded from above by 



J2([AA^]^<Sh'N^ (69) 



as can be shown by using a symmetrization procedure with subsequent ap- 
phcation of the spectral theorem, similar to the procedure outlined in the 
appendix of |jlT| . 



For the Hamiltonian-dependent double commutator, one finds as an upper 
bound 

[IC H]_ ,C^]_)< ANdh^ {\BoM{T, Bo)\ + 2gP) , (70) 



where the constant q is taken to reflect the fact that the hopping constants 



tf fulfill the convergence criterion 



A^LL^n. 1 = ^<oo. (71) 

7/3 nk 

Inserting these results into the (fc-summed) Bogoliubov inequality, and 
making the transition to the thermodynamic limit as in (^Ol), one arrives at 
an inequality for the layer magnetization which is similar, though not 
identical with, the one in (^21): 



M',{T, Bo) < i—. (72) 

V \BoM{T,Bo)\ ) 

the main difference being the factor d (=number of layers) on the right- 
hand side. Interestingly, the upper bound is proportional to the inverse 
temperature 13 and the number of layers d. While for any finite (3 (T 7^ 0), 
and for any finite number of layers {d < 00) a phase transition is ruled out 
because of the divergence of the denominator as i?o — > 0, the possibility of 
a phase transition opens up if either j3 01 d are infinite, i.e. when a two- 
dimensional system is considered at T = 0, or when the system becomes 
truly three-dimensional (i.e. infinitely extended in all dimensions). 

As an example for a superconducting phase transition, we shall now 
consider a pairing mechanism for Hubbard films. Following the discussion in 
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PB| , we restrict our attention to the order parameter 



= ]^ EE ^-"■■"•(4,4,, 

a i 

^ \ll^o.{K) (73) 

i 

which measures the breakdown of U(l) symmetry due to local on-site pair- 
ing {'^a{,K) is again the layer-dependent equivalent to the bulk "^{K)). The 
ordering wave vector allows one to distinguish between different types of 
pairing: for K = {] one would have s-wave pairing, for K ^ (generalized) 
//-pairing. The Hubbard Hamiltonian with an appropriate U(l) symmetry- 
breaking contribution of order A is then given by 



ija/3cr iaa 



-^E(^^^«^"^''''^'+^r.e+^^-^^') (74) 

ia 

where we have introduced the operators 

Via = 4a1'4al C^'^) 
Via = CiaiCia^ (76) 

via = 2 (^^°T + - 2) • (77) 
With the operators A and C chosen as 

C = 5^<(^) (78) 

7 

A = r^^{-k-K), (79) 
the calculation of the (anti-)commutators for the Bogoliubov inequality re- 
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suits in 



EE' 



,i(/c+i^)-(ij,-Rj) / r„+ 



i 

[[C,H]^,C^_) < Ndqk^ + 2XNd ^{K) . 



(81) 
(82) 



Thus, in this case as well, one gets a result of the form 



(2) 
271 



< 



(3d 



In ( 1 + 



(83) 



which excludes a finite value of the layer-specific quantity '^a, again due to 
the divergence of the denominator in the limit A — *^ 0. Thus, no pairing 
transition of the proposed kind and hence no corresponding superconducting 
phase transition can occur in Hubbard films, provided the number of layers 
is finite and the temperature is non-zero. 

The dependence on d of the upper bound in the film cases, confirms 
a statement by Fisher which is based on his work with Jasnow |5D|, 
where the upper bound scales with the extension of the "box" A in which 
the physical system is embedded. Note, however, that the two approaches 
start from opposite ends: Fisher and Jasnow impose external geometrical 
constraints, whereas the approach presented here (and the related approach 
of ref. [|lT|) includes the film geometry in the Hamiltonian from the very 
beginning. 



3.4 Fractal lattices 

The discussion in the previous section points to a deeper connection between 
the geometry of the sample and the interaction between its constituents. 
Hattori et al. ||4^ have pointed to an interesting aspect of this relation, and 
we shall briefly summarize their discussion. 
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In the theory of phase transitions, and particularly in the versions of the 
Mermin- Wagner theorem discussed so far, one typically deals with a spin 
system on a translationally invariant lattice which is embedded into a Eu- 
clidean space of integer dimension {D = 1,2, or 3). It turns out that the 
conditions for the (non-)existence of a phase transition are governed by the 
dimension of this Euclidean space. More generally, the critical properties 
and the scaling limits of such a many-particle system crucially depend on 
the Euclidean metric associated with the dimensionality: Correlation func- 
tions of the form are Euclidean invariant functions in the scaling 
limit. This, however, poses a question which Hattori et al. put forward in 
the following way: "Why do the spins on the lattice 'know' this natural em- 
bedding into Euclidean space that should govern their critical phenomena?" 
After all, the only geometrical structure that spins can feel is the one given 
by the Hamiltonian, i.e. by the "network structure of the interaction (the 
kinetic term, in terms of field theory)." |^ 

As long as one is concerned with highly regular, i.e. translationally in- 
variant lattices, this duality between the structure of the interaction (as 
contained in the Hamiltonian) and Euclidean dimension remains somewhat 
hidden. One instance where it surfaces, however, is the complementarity, 
mentioned in the last paragraph of the previous section, between Fisher and 
Jasnow's approach (which suggests that the upper bound in the Mermin- 
Wagner approach is ~ A, i.e. proportional to the linear extension of the 
embedding box") and our calculation which proves that the upper bound in 
a film system is indeed proportional to the number of layers d (by virtue of 
the structure of the Hamiltonian). 

The question becomes highly relevant for irregular lattices, where many of 
the techniques (e.g. Fourier transforms, or notions of "momentum" vectors) 
cannot be applied any longer. As a step towards discussing such irregular 
systems, Hattori et al. develop a Gaussian field theory on a general network; 
the authors succeed in giving a definition of the spectral dimension of general 
networks in terms of the critical behaviour of a spin system, and show how, 
for a restricted class of networks, the spectral dimension can be determined 
in practice. While a detailed discussion of their results would be far beyond 
the scope of this paper, an application by Cassi |^ is of interest. By ap- 
plying Hattori et al.'s results to the case of fractal and disordered lattices 
whose spectral dimension is less or equal to 2, it is shown that at finite tem- 
perature no spontaneous magnetization can exist for the classical 0{n) and 
the ferromagnetic (quantum) Heisenberg model. This is due to the infrared 
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divergencies displayed by the corresponding Gaussian models and confirms 
once again the wide applicability of the Mermin- Wagner theorem. 



4 On the validity of the Mermin- Wagner the- 
orem in approximative theories 

Surprisingly few authors have commented on the relevance of the Mermin- 
Wagner theorem for methods of approximation, some of which consistently 
predict finite-temperature phase transitions in one or two dimensions, even 
when applied to many-body models that obey the Mermin- Wagner theorem. 
Sometimes the breaking of the Mermin- Wagner theorem has been taken to 
be an advantage; Sukiennicki and Wojtczak, for example, believe "that the 
molecular field approximation of the Hubbard model approaches the phys- 
ical reality (magnetic order in thin films does really exist) better than the 
Hubbard model in its exact form." |Q While this may seem prima facie rea- 
sonable, it is difficult to maintain this view in the light of the advances that 
have been achieved within the (unmodified) Hubbard model by disposing 
with simple molecular field approximations and adopting instead more re- 
fined methods of approximation. If one aspires to formulate criteria for the 
reliability of one's methods of approximation (other than justifying them 
in an ad hoc fashion), one will need to shed some light on the relation 
between the underlying many-body model, the exact results known about 
it, and the mechanism of approximation. It would seem unsatisfactory to 
rely blindly on the method of approximation introducing exactly the kind 
of symmetry-breaking needed to reflect physical reality. If it turns out that 
certain classes of approximation methods consistently produce reasonable 
results for quasi-2D cases, this needs explanation. One could even turn the 
argument around and question methods of approximation for the 3D case in 
case they violate the Mermin- Wagner result in two dimensions. As Walker 
puts it: "[Ajpproximation schemes that have been applied to real solids can 
equally well be applied to one- and two-dimensional solids. If these approxi- 
mation schemes predict the occurrence of spontaneous magnetization in one 
and two dimensions as well as in three dimensions [...], the validity of these 



predictions in three dimensions should clearly be investigated more fully." |22 
Considerations of this sort set a possible agenda for future research into 
the relevance of exact results for the application of approximation schemes. 
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In the rest of this section, we shall present a fairly elementary discussion of 
the Heisenberg model within different methods of approximation. In partic- 
ular, we shall contrast two possible ways of testing "Mermin- Wagner type" 
behaviour. 

Consider the spin-1/2 anisotropic Heisenberg model given by 

^ = - E iJ^^ (^^^i + ^'^i ) + 4^^^!) (84) 

n.n.ij 

where = eJij > Jij and the summation index ^'"n.n." refers to summation 
over nearest neighbours only. This kind of anisotropy in spin-space favouring 
uniaxial ordering is widely used and should, even in 2D, break the symmetry 
of the corresponding isotropic model. 

It has sometimes been claimed that a modification of mean-field 



theory, known as Onsager reaction-field theory, can reproduce the Mermin- 
Wagner theorem in the case of the Heisenberg model. According to the On- 
sager theory the orienting part of the mean field acting on a given spin must 
not include that part of the contribution from the spins in the vicinity which 
is due to their correlation with the given spin under consideration. Thus, 
one must stipulate that the full mean field decomposes into two independent 
contributions, the correlated and the uncorrelated one, called reaction field 
and cavity field, respectively. (For a detailed discussion see ^^.) Singh ||5^ 



has succeeded in applying Onsager's reaction-field idea to the anisotropic 



Heisenberg model as given by (|8^). It turns out that the Onsager reaction- 
field result for the critical temperature, T^"-^^ is related to that of ordinary 
mean-field (MF) theory by the equation 



T 

c 



Ons 



TMF ^ 



1 



^5) 



k V J'-'iK), 



where K denotes the ordering wave vector and J^{k) is the standard Fourier 
transform of the exchange integrals Jij. From (|85|) , however, we see that 
while Onsager reaction-field theory does modify the ordinary mean-field re- 
sult, it is independent of the anisotropy assumed above by setting J^^ = eJij, 
since the anisotropy parameter e will drop out of the fraction J^{k) / J^{K) 
in the denominator. Effects of dimensionality that stem from the summation 
are, on the other hand, retained, such as the denominator's logarithmic 
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divergence in two dimensions, which rules out a phase transition. This com- 
phance with the Mermin- Wagner theorem, however, is rendered spurious by 
the fact that no anisotropy effects are captured by the Onsager approach. 

It is worthwhile to explore whether more sophisticated methods of ap- 
proximation are more successful in accounting for effects due to (reduced) 
dimensionality and (possibly anisotropic) interaction. For the Tyablikov pro- 
cedure of decoupling higher-order Green's functions within the equation of 



motion scheme, this is indeed the case. As has been shown in [^, the Tyab- 
likov procedure accords with the Mermin- Wagner theorem in two dimensions: 
Only when an external field is applied, can there be a finite value of the mag- 
netization at non-zero temperature. The proof is based on a series expansion 
of the (implicit) equation for the expectation value (S^), 



hs 




where the energy E{k) is given by E{k) = 2h (5^) ^Jq — J{k)j +gjfiBBQ. As a 

result of the Tyablikov approximation, (|86|) is taken to be uniform for all spins 
in the lattice, which seems to limit the discussion to the ferromagnetic case. 
However, it has also been shown that an analogous result holds for the case of 
spontaneous sublattice magnetization in ABAB-type antiferromagnets. The 
corresponding formula for the Neel temperature in that case is of the following 
form: 

T^' = ' , • (87) 



1- 



It is straightforward to show that in the thermodynamic limit the sum in the 
denominator diverges in two dimensions while remaining finite in three di- 
mensions. What remains to be shown is that the Tyablikov method can make 
sense of anisotropy effects, i.e. perform better than the Onsager reaction-field 
theory in accounting for ordering due to anisotropy. One such test is whether 
or not the method can correctly predict the Mermin- Wagner theorem if one 
starts from the three-dimensional case and gradually turns off interlayer cou- 
pling. If we choose the interaction parameters Jij so that (for a fixed lattice 
site i) Jij = J|| if sites i and j lie within a plane, and Jij = J± if j is located 
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in the direction orthogonal to this plane, then, by varying the anisotropy 
parameter 



e := ^ (88) 

we can change the coupling continuously from quasi two-dimensional {e = 
0) to (isotropic) three-dimensional {e = 1). One can then, indeed, prove 
that in the quasi-2D limit of vanishing interlayer coupling {e — ^ 0) the Neel 
temperature (as calculated from the Tyablikov method) vanishes as 

(89) 



|ln(e)|' 

in accordance with the Mermin- Wagner theorem. In its original form, the 
Tyablikov method applies to the S = 1/2 case only, and so, by definition, do 
the results derived so far. An extension of the results to higher spins can. 



however, be given in a straightforward way |^ . 

In this section, we have contrasted two approximation schemes, the On- 
sager reaction-field theory and the Tyablikov method. The failure of reaction- 
field theory to capture the effect of anisotropic coupling suggests that it does 
not adequately capture the behaviour of the Heisenberg model as far as the 
existence, or absence, of a phase transition is concerned. The Tyablikov 
method, despite its still rather simple random-phase characteristics, seems 
to fare better: It can be analytically shown to produce all the relevant limit- 
ing cases of possible "Mermin- Wagner type behaviour," and it also correctly 
reproduces the qualitative change associated with the transition from three 
to two dimensions by introducing anisotropic (and eventually vanishing) in- 
terlayer coupling. 



5 Summary 

In this paper, we have presented a survey of recent and established results 
concerning the application of Bogoliubov's inequality to the theory of phase 
transitions. This includes the classic papers by Hohenberg, Mermin, and 
Wagner as well as a host of other proofs concerning a variety of different 
many-body models, order parameters and system geometries. In particular, 
new proofs for the absence of superconducting long-range order in Hubbard 
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films, and of magnetic long-range order in Periodic Anderson Model films 

have been presented. The complementarity between geometric constraints 
on the one hand and anisotropy on the other hand, has been sketched in 
theoretical terms, and has been applied to two specific methods of approxi- 
mation, Onsager reaction-field theory and the Tyabhkov method. The two 
methods discussed here, can only give an indication as to which effects and 
limiting cases should be considered, if one attempts to characterize a method 
in terms of its accordance with the Mermin- Wagner theorem. The more gen- 
eral problem of the relevance of exact results for the numerical or approx- 
imative treatment of many-body models, as characterized in section 4, will 
almost certainly remain and in our view deserves closer attention. 
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